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Abstract 

We report results from the Eddy Experiment, where a synchronous GPS receiver pair was flown on an 
aircraft to collect sampled LI signals and their reflections from the sea surface to investigate the altimetric 
accuracy of GNSS-R. During the experiment, surface wind speed (U10) was of the order of 10 m/s, and 
significant wave heights of up to 2 m, as discussed further in a companion paper. After software tracking 
of the two signals through despreading of the GPS codes, a parametric waveform model containing the 
description of the sea surface conditions has been used to fit the waveforms (retracking) and estimate the 
temporal lapse between the direct GPS signals and their reflections. The estimated lapses have then been 
used to estimate the sea surface height (SSH) along the aircraft track using a differential geometric model. 
As expected, the precision of GNSS-R ranges was of 3 m after 1 second integration. More importantly, the 
accuracy of the GNSS-R altimetric solution with respect to Jason- 1 SSH and in situ GPS buoy measure- 
ments was of 10 cm, which was the target with the used experimental setup. This new result confirms the 
potential of GNSS-R for mesoscale altimetric monitoring of the ocean, and provides an important milestone 
on the road to a space mission. 
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1. Introduction 

Several Global Navigation Satellite System 
(GNSS) constellations and augmentation systems 
are presently operational or under development, 
such as the pioneering US Global Positioning Sys- 
tem (GPS), the Russian Global Navigation Satellite 
System (GLONASS) and the European EGNOS. In 
the next few years, the European Satellite Naviga- 
tion System (Galileo) will be deployed, and GPS 
will be upgraded with more frequencies and civil- 
ian codes. By the time Galileo becomes opera- 
tional in 2008, more than 50 GNSS satellites will 
be emitting very precise L-band spread spectrum 
signals, and will remain in operation for at least a 
few decades. Although originally meant for (mil- 
itary) localization, these signals will no doubt be 
used within GCOS/GOOS* in many ways (e.g., at- 
mospheric sounding). We focus here on the bud- 
ding field known as GNSS Reflections, which aims 
at providing tools for remote sensing of the ocean 



surface (especially sea surface height and roughness) 
and the atmosphere over the oceans. 

This paper reports a new development of the 
GNSS-R altimetric concept (PARIS). The PARIS 
concept (Passive Reflectometry Interferometric Sys- 
tem 1121 1 addresses the exploitation of reflected 
Global Navigation Satellite Systems signals for al- 
timetry over the oceans. Ocean altimetry, the mea- 
surement of Sea Surface Height (SSH), is indeed 
one of the main applications of the GNSS-R pas- 
sive radar concept. GNSS-R can give automatically 
integrated measurements in the GNSS reference sys- 
tem. In addition, this technique can provide the un- 
precedented spatio-temporal samplings needed for 
mesoscale monitoring of ocean circulation. It is at 
mesoscale where phenomena such as eddies play 
a fundamental role in the transport of energy and 
momentum, yet current systems are unable to probe 
them. 



* Global Climate Observing System/Global Ocean Observing System. 
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Many GNSS-R altimetry and scatterometry ex- 
periments have been carried out to date, and the 
list continues to grow thanks to dedicated efforts 
in Europe and the US. GNSS-R experimental data 
has now been gathered from Earth fixed receivers 
(lUIISlEOlEl among others), aircraft (|9| HE EH HQ 
among others), stratospheric balloons (|6, 4| among 
others), and from space platforms (| 10 1 among oth- 
ers). This experimental work is converging to a uni- 
fied understanding of the GNSS-R error budget, but 
so far these experiments have focused on waveform 
modeling and short term ranging precision. None to 
date have attempted to retrieve a mesoscale altimet- 
ric profile as provided by monostatic radar altimeters 
such as Jason- 1. 

In the four main sections of this paper we re- 
port PARIS altimetric processing results using data 
from the 09-27-2002 Eddy Experiment, carried out 
in the frame of the European Space Agency "PARIS 
Gamma" contract. The first section addresses the 
issue of tracking the direct and reflected GPS sig- 
nals, which consist in appropriately placing the de- 
lay and Doppler gating windows and in despread- 
ing the GPS signals by means of correlation with 
clean replicas. Tracking produces incoherently aver- 
aged waveforms (typically with a cadence of 1 sec- 
ond). The extraction of the information needed for 
the altimetric algorithm from the waveforms is de- 
scribed in the second section. This is the retracking 
step, and it yields to the so-called measured tempo- 
ral lapse (or lapse, for short) between the direct and 
reflected signal. In the third section, the altimetric 
algorithm (producing the Sea Surface Height esti- 
mates) is described and, finally, results are presented 
in the fourth section. 

2. Data collection and pre-processing 
2.1. Data set 

The GNSS-R data set was gathered during an air- 
borne campaign carried out by Starlab in Septem- 
ber 2002. The GPS/INS (Inertial Navigation Sys- 
tem) equipped aircraft overflew the Mediterranean 
Sea, off the coast of Catalonia (Spain), northwards 
from the city of Barcelona for about 150 km (Fig- 
ureQ. This area was chosen because it is crossed by 
a ground track of the Jason- 1 altimeter (track num- 
ber 187). The aircraft overflew this track during the 
Jason- 1 overpass, for precise comparison. In addi- 

tThe buoy campaign was carried out by the Institut d'Estudis Es 
* STARLIGHT, also described in 0. 



tion, a GPS buoy measurement of the SSH on a point 
along the same track was obtained, basically to vali- 
date the Jason- 1 measurement^. 

During the experiment, the aircraft overflew 
twice the Jason- 1 ground track, both in the South- 
North direction and back, gathering about 2 hours 
of raw GPS-L1 IF data sampled at 20.456 MHz (see 
11171 for the experimental setup). In this paper, we 
focus on the data processing of the ascending track, 
i.e., from PI to P2 in Figure^] with the Jason-1 over- 
pass happening roughly in the middle of the track. 

2.2. Tracking GNSS-R signals 

Altimetry with GNSS-R signals is based on the mea- 
surement of the temporal lapse between the time of 
arrival of the direct GNSS signal and the time of ar- 
rival of the same signal after its reflection from the 
target surface. Successful tracking of both signals is 
the first step for altimetric processing. 

Under general sea conditions, GPS signals re- 
flected from a rough sea surface cannot be tracked by 
a standard receiver, because of the signal corruption 
due to the reflection process itself (in GPS terminol- 
ogy, the signal is affected by severe multipath from 
the sea clutter). For this reason, a dedicated soft- 
ware receiver has been developed^. The high level 
block diagram of this receiver is shown in Figure|2] 
The processor is composed of two sub-units, one 
for each signal. The unit which processes the direct 
signal — the master unit — uses standard algorithms 
to track the correlation peak of the signal, both in 
time and frequency. The unit which processes the 
reflected signal — the slave unit — performs correla- 
tions blindly, with time delay and frequency settings 
which depend on those of the master unit. 

One of the most relevant tracking parameters is 
the temporal span of the correlations, i.e., the coher- 
ent integration time used to despread the GPS signal. 
The coherent integration time was set here to 10 mil- 
liseconds: it was verified that with this value the ra- 
tio of the correlation peak height and the out-of-the- 
peak fluctuations achieved a maximum. In practical 
terms, an integration time of 10 milliseconds sim- 
plifies the tracking process, as the duration of this 
time interval is a sub-multiple of the GPS navigation 
bit duration (with the 50 Hz navigation code in LI). 
Moreover, we believe that longer correlation times 
provide some protection from aircraft multipath by 

icials de Catalunya (IEEC). 
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Figure 1: Flight trajectory. 



mimicking a higher gain antenna (a belief supported 
by tests with shorter integration times). 

3. Retracking the waveforms 

Once a correlation waveform is obtained for both 
the direct and reflected signals, the lapses can be es- 
timated. We emphasize that this is not as trivial as 
considering the maximum sample of each waveform 
or the waveform centroid, for instance, as the bistatic 
reflection process deforms severely the signals and, 
in general, such distortions will displace the wave- 
form maximum or centroid location. Moreover, the 
sampling rate of 20.456 MHz, while much higher 
than the Nyquist rate for the C/A code, is equivalent 
to an inter-sample spacing of 49 ns — about 15 light- 
meters. This coarseness in the lapse estimation is 
not sufficient to reach the final altimetric precision 
target. The main challenge for GNSS-R using the 
GPS Coarse Acquisition (C/A) code is to provide 
sub-decimetric altimetry using a 300 m equivalent 
pulse length, something that can only be achieved 
by intense averaging with due care of systematic ef- 
fects. For reference, pulse lengths in monostatic al- 
timeters such as Jason- 1 are almost three orders of 
magnitude shorter. 



For these reasons, the temporal position of each 
waveform (the delay parameter) is estimated via a 
Least Mean Squares model fitting procedure. This is 
the so-called retracking process, an otherwise well- 
known concept in the monostatic altimetric com- 
munity. The implementation of accurate waveform 
models (for direct and reflected signals) is funda- 
mental to retracking. Conceptually, a retracking 
waveform model allows for the transformation of the 
reflected waveform to an equivalent direct one (or 
vice-versa), and a coherent and meaningful compar- 
ison of direct and reflected waveforms for the lapse 
estimation. 

3.1. Waveform model 

The natural model for retracking of the direct signal 
waveforms is the mean autocorrelation of the GPS 
C/A code in presence of additive Gaussian noise, 
which accounts mainly for the receiver noise. As 
far as the reflected signal is concerned, the model is 
not so straightforward. In fact, the reflection process 
induces modifications on the GNSS signals which 
depend on sea surface conditions (directional rough- 
ness), receiver-emitter-surface kinematics and ge- 
ometry, and antenna pattern. Among all these quan- 
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Figure 2: Simplified block diagram of the GNSS-R tracking concept. 



tities, the least known ones are those related to the 
sea surface conditions. In principle, these quantities 
should be considered as free parameters in the model 
for the reflected signal waveform and estimated dur- 
ing the retracking process along with the delay pa- 
rameter of the waveform. 

On the basis of the two most quoted models in 
the literature for bistatic sea-surface scattering (| 14 1 
and 1 21 1), we have developed an upgraded waveform 
model for the reflected signal. This new model, as 
the two aforementioned, is based on the Geometric 
Optics approximation of the Kirchhoff theory — that 
is to say, with a tangent plane, high frequency and 
large elevation assumption, which is reasonable for 
the quasi-specular sea-surface scattering at L-band 
(' 1171 '). The main characteristics of this model are 
the following: 

• a fully bistatic geometry (as in [21 1, but not in 

O), 

• the description of the sea surface through 
the L-band Directional Mean Square Sloped 
(DMSSl) (as in ED), and 

• the use of a fourth parameter, the significant 
wave height (SWH), to describe the sea sur- 
face (as in 1 14 1, but not in 1211 ). 

We have checked that the impact of SWH mismod- 
eling in our case is negligible, since the GPS C/A 
code equivalent pulse width is about 300 meters. 
Nonetheless, we emphasize that all potential sources 

§ See 1 8 for a discussion on the role of wavelength in the definiti 



of systematic effects must be considered. We fore- 
see a higher and non-negligible impact of SWH if 
the GPS P-code (the Precision code) or similar codes 
in Galileo are used, since they are an order of mag- 
nitude shorter. 

3.2. Inversion scheme 

The retracking process has been performed through 
a Least Mean Squares fit of the waveforms with 
the models described. Because of the speckle 
noise affecting the reflected signal waveforms, the 
fit has not been performed on each complex wave- 
form obtained from the 10 ms correlations. Rather, 
these waveforms have first been incoherently aver- 
aged: the average of the magnitude of a hundred 
10 ms waveforms has then been used to perform the 
inversion — i.e., 1 second incoherently averaged real 
waveforms have been generated for retracking. In 
this way, reflected/direct temporal lapses have been 
produced at 1 Hz rate. 

In both cases, the fit of the waveform has been 
performed over three parameters: the time delay, 
a scaling factor and the out-of-the-peak correlation 
mean amplitude. 

The geophysical parameters that enter in the 
model of the reflected signal waveform have not 
been jointly estimated here. These parameters have 
been set to some reasonable a priori value obtained 
from other sources of information (Jason- 1 for wind 
speed, ECMWF for wind direction) or from theory 
(for the sea slope PDF isotropy coefficient). For con- 

tofDMSS. 
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venience, we describe the sea surface state using a 
wind speed parameter, wind direction and the sea 
slope PDF isotropy coefficient. Using a spectrum 
model, these parameters can be univocally related 
to DMSSl with the assumption of a mature, wind- 
driven sea (the sea spectrum in 1 5 1 has been used in 
this case). We emphasize that DMSSl is the actual 
parameter set needed in the reflection model under 
the Geometric Optics approximation. 

Concerning the reflected signal waveform, re- 
tracking has been performed using only the lead- 
ing edge and a small part of the trailing edge, since 
the trailing edge is more sensitive to errors in the 
input parameters (including geophysical parameters 
and antenna pattern). 

4. The altimetric algorithm 

The output of the retracking process is the time 
series of measured lapses. The next step is finally 
SSH estimation. In order to solve for this quantity, 
we have used a differential algorithm: the classi- 
cal PARIS equation (see 1 12 1) has not been directly 
used. Instead, a model for the lapse over a refer- 
ence surface near the local geoid has been built, and 
the difference of this reference lapse and the mea- 
sured one has been modeled as a function of the 
height over the reference surface. We call this the 
differential PARIS equation (Equation [0- We note 
that the aircraft INS has been used to take into ac- 
count the direct-reflected antenna baseline motion 
and that we have also included both dry and wet tro- 
pospheric delays in the model by using exponential 
models for them with different scale heights and sur- 
face values derived from surface pressure measure- 
ments and surface tropospheric delays obtained from 
ground GPS and SSM/I^). 
The differential PARIS equation writes 

A DM = Ac — Am = 2 Sh sin(e) + b, (1) 

where 

• Arj is the measured lapse, in meters, as esti- 
mated from the data, 

• Am is the modeled lapse, in meters, based on 
an ellipsoidal model of the Earth, GPS con- 
stellation precise ephemeris, aircraft GPS/INS 
precise kinematic processing, and a tropo- 
spheric model, 

^Special Sensor Microwave Imager, a passive microwave radio! 
Program. 



• 5h is the normal offset between the sea surface 
and the (model) ellipsoid surface, 



• e is the GPS satellite elevation angle at the 
specular point of reflection, over the ellipsoid, 
and 



• b is the hardware system bias. 

The precision obtained after 1 -second of incoher- 
ent averaging in the estimation of Aom using this 
approach is displayed in Table 1. For each PRN 
number, the root mean squared error of the 1 -second 
lapse is shown (in meters). It is roughly of 3 m. This 
noise level is as expected from the C/A code bistatic 
ranging in our experimental setup (antenna gain, al- 
titude, etc.) and consistent with earlier experiments. 



5. Results 

The algorithm outlined in the previous section 
has been used with data from the three better be- 
haved GPS satellites. The other two visible satel- 
lites appeared to be severely affected by aircraft- 
induced multipath (probably due to the landing gear 
and wing surfaces). A SSH profile has been retrieved 
for each satellite and the average value of the three 
profiles is shown in Figure [5] along with the Jason- 
1 SSH, Mean Sea Level (MSL) and one SSH value 
obtained with the control buoy. This solution has 
been obtained setting a model wind speed of 10 m/s 
(values provided by Jason- 1 along the track vary be- 
tween 9 and 13 m/s), wind direction of degrees 
North (values provided by ECMWF vary between 
30 and -20 deg North), and sea slope PDF isotropy 
coefficient equal to 0.65 (the theoretical value for a 
mature, wind-driven sea according to 0). It is im- 
portant to underline that the use of constant values 
for the geophysical parameters along the whole track 
(more than 120 km) induces non-linear errors on the 
final altimetric estimation. Nonetheless, the bias of 
the final solution with respect to the SSH (the error 
mean) is 1 .9 cm while the root mean error is 10.5 cm. 

:ter flown aboard United States Defense Meteorological Satellite 
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PRN 


Complete track 


Beginning of the track 


8 


3.5 m 


2.7 m 


24 


3.4 m 


2.8 m 


27 


2.9 m 


2.7 m 



Table 1: Precision in the estimation of the time lapses (root mean squared error of the lapses, in meters). 
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Figure 3: Final altimetric result. The red line represents the average altimetric result of the three GPS satellite analyzed filtered over a 
20 km — i.e., 400 seconds, the aircraft speed being about 50 m/s. The black line represents the Jason-1 SSH, while the green dashed 
line represents the MSL. The blue diamond is the SSH measurement obtained from the reference buoy. 



6. Conclusion 

The Eddy Experiment has validated the use of 
PARIS as a tool for airborne sea surface height mea- 
surements, providing both the precision and accu- 
racy predicted by earlier experimental and theoreti- 
cal work. We have observed that the use of a wave- 
form model for the reflected signal, based on geo- 
physical parameters describing the sea surface con- 
ditions, is essential for the accuracy of the altimetric 
solution — a fact which may explain earlier results 
in which no geophysical retracking was performed 
(e.g., [15 1). The accuracy achieved by our algorithm 
is of the order of 1 decimeter, but we expect that fur- 
ther analysis and refinements, such as the inclusion 
of DMSSl parameters in the inversion procedure, 
will improve these numbers. 

Our sensitivity analysis has also shown that the 
altitude of this flight was not optimal for GNSS-R 
altimetry and made the experiment more prone to 



aircraft multipath problems. A higher altitude flight 
would lead to a smaller angular span of the reflect- 
ing area on the sea surface, thus reducing the impact 
of geophysical parameters, antenna pattern and air- 
craft multipath on the retracking process of the lead- 
ing edge, making the overall altimetric solution more 
robust. Higher altitudes are also needed to better un- 
derstand the space-borne scenario. 

We would like to emphasize that GNSS-R sig- 
nals can be profitably used also for scatterometric 
measurements (i.e., speculometry, from the Latin 
word for mirror, "speculo", see |16|). In a par- 
allel paper (|8|), the inversion of GNSS-R Delay- 
DopplerMaps for sea-surface DMSSl estimation is 
presented for the same data set. The next step is to 
merge the altimetric and speculometric processing in 
an attempt to provide an autonomous GNSS-R com- 
plete description of the sea — topography and surface 
conditions. 

We believe the Eddy Experiment is an important 
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milestone on the road to a space mission. We un- 
derline that the obtained precision and accuracy are 
in line with earlier experiments and theoretical error 
budgets (see, e.g., 1 10 1). We note that the same error 
budgets have been used to investigate and confirm 
the strong impact of space-borne GNSS-R altimetric 
mission data on mesoscale ocean circulation mod- 
els ( 11811191 1. Further analysis of existing datasets 
(which could be organized in a coordinated database 
for the benefit of the community) and future experi- 
ments at higher altitudes will continue to refine our 
understanding of the potential of this technique. 
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